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A status report is made of an on-going full QCD study on the CP-PACS aiming at a comparative analysis of 
the effects of improving gauge and quark actions on hadronic quantities and static quark potential. Simulations 
are made for four action combinations, the plaquette or an RG-improved action for gluons and the Wilson or SW- 
clover action for quarks, at a" 1 « l.l-1.3GeV and m n /m p w 0.7-0.9. Results demonstrate clearly that the clover 
term markedly reduces discretization errors for hadron spectrum, while adding six-link terms to the plaquette 
action leads to much better rotational symmetry in the potential. These results extend experience with quenched 
simulations to full QCD. 



1. Introduction 

With the progress in recent years of quenched 
simulations of QCD, deviations of the quenched 
hadron spectrum from experiment are being un- 
covered. For heavy quark systems precise calcu- 
lations with NRQCD have shown that the fine 
structure of quarkonium spectra can be repro- 
duced only if sea quark effects are taken into ac- 
count Q. For the light hadron sector several re- 
ports have been made that strange quark mass 
cannot be set consistently from pseudo scalar and 
vector meson channels in quenched QCD ^ ||. 
Most recently results of an extensive quenched 
simulation on the CP-PACS indicate that there 
is a systematic disagreement in the spectrum of 
baryons || . Clearly the time has come to bolster 
efforts toward full QCD simulations. We have re- 
cently started an attempt in this direction using 

"Talk presented by K. Kanaya at the International Work- 
shop on "LATTICE QCD ON PARALLEL COMPUT- 
ERS", 10-15 March 1997, Center for Computational 
Physics, University of Tsukuba. 



the CP-PACS computer §. 

Full QCD simulations are, however, extremely 
computer time consuming compared to those of 
quenched QCD. Even with the TFLOPS-class 
computers that are becoming available, high 
statistics studies, indispensable for reliable re- 
sults, will be difficult for lattice sizes exceeding 
32 3 x 64. Since a physical lattice size of L w 2.5- 
3.0fm is needed to avoid finite-size effects 
the smallest lattice spacing one can reasonably 
reach will be a -1 » 2GeV. Hence lattice dis- 
cretization errors have to be controlled with sim- 
ulations carried out at an inverse lattice spac- 
ing smaller than this value. This will be a diffi- 
cult task with the standard plaquette and Wilson 
quark actions since discretization errors are of or- 
der 20-30% even at a" 1 w 2GcV Q. This leads 
us to consider improved actions for our simulation 
of full QCD. 

Studies of improved actions have been widely 
pursued in the last few years. Detailed tests of 
improvement for hadron spectrum, however, have 
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Table 1 

Simulation parameters of full QCD runs on a 
12 3 x 32 lattice. Numbers in parenthesis denote 
number of configurations stored for static quark 
potential analyses. 
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been mostly carried out within quenched QCD 
(see, e.g., Refs. fl2]-|l5"|), and only a few studies 
are available for the case of full QCD ju|. In 
particular a systematic investigation of how vari- 
ous terms added to the gauge and quark actions, 
taken separately, affect light hadron observables 
has not been carried out in full QCD. We have 
decided to undertake such a study as the first 
subject of our full QCD program. In this arti- 
cle we report preliminary results of this on-going 
attempt. 

2. Choice of action and simulation param- 
eters 

Improving the standard plaquette action for 
gluons requires the addition of Wilson loops of six 
links or more in length. The precise forms of the 
added terms and their coefficients differ depend- 



ing on the principle one follows for improvement. 
In our study we choose an action given by 

5 f =(K c °^^ ixi+ci ^^ ix2 )' (i) 

with cO = 1 — 8ci and c\ = —0.331, which was ob- 
tained by a renormalization group treatment fTi| . 
The quenched static quark potential calculated 
with this action exhibits good rotational symme- 
try and scaling already at a -1 » lGeV|ll|, sim- 
ilar to those observed for tadpole-improved and 
fixed point actions pjj , p0| . 

For improving the quark action we take the 
clover improvement due to Sheikholeslami and 
Wohlert @ defined by 

Dg = D%+ 5 xy c sw K l^uF^ (2) 
with the standard Wilson matrix D^f y given by 

+ (1 + 7mRV<Wa}= ( 3 ) 

For the clover coefficient cswi we adopt the 
meanfield improved value p2| ]: c$w = -P -3 / 4 
with P the plaquette average. 

We carry out a comparative study of the light 
hadron spectrum and static quark potential for 
four action combinations, choosing either the pla- 
quette (P) or the rectangular action above (R) for 
gluons and either the Wilson (W) or clover action 
(C) for quarks. We expect the extent of improve- 
ment to be clearer at a coarser lattice spacing. We 
therefore attempt to tunc the coupling constant 
so that the lattice spacing equals a -1 ~ lGeV. 

Our simulations are carried out for two flavors 
of quarks, mostly on a 12 3 x 32 lattice with addi- 
tional runs on an 8 3 x 16 lattice to estimate pa- 
rameters, including a self-consistent value of csw- 
We employ the hybrid Monte Carlo algorithm to 
generate full QCD configurations at two or three 
values of K corresponding to m^/m p w 0.7-0.9. 
The molecular dynamics step size is chosen to 
yield an acceptance of 80-90%. After thermaliz- 
ing for 100-200 trajectories we generate 1000-1500 
trajectories. 
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Figure 1. (a) mjq jm p and (b) m&/m p as a function of (m^/mp) 2 for various combinations of the action. 
Stars in (b) are experimental points corresponding to A(1232)/p(770) and Q(1672)/</>(1020). 



We measure hadron propagators every 5 trajec- 
tories, using point and smeared sources and point 
sinks following the method of our quenched study 
H). The static quark potential is calculated on 
a subset of configurations used for hadron prop- 
agator measurement. The smearing technique 
of Ref. |23| is employed choosing the number of 
smearing steps and fitting ranges from experience 
in Ref. |jl8|| . Errors are estimated by a single- 
elimination jackknife procedure. Simulation pa- 
rameters of our runs and the number of configu- 
rations used for the spectrum and potential mea- 
surements are summarized in Table [l| 

3. Light hadron masses 

Our main results for the effect of improved ac- 
tions on hadron masses are displayed in Fig. [I] 
in which the ratio m^r/mp and m&/m p are plot- 
ted as a function of (m^/mp) 2 for the four ac- 
tion combinations. The solid curves represent 
the well-known phenomenological mass formula 
(24). The inverse lattice spacing estimated from 
the p meson mass in the chiral limit is in the range 
a- 1 ps l.l-1.3GeV (see Table |). 

For the standard action combination P-W, the 
ratios are well above the phenomenological curve 
as may be expected at such a large lattice spac- 
ing. When we improve the gauge action (the R-W 
case), the data points come closer to the curve. 
By far the most conspicuous change, however, is 



Table 2 

Results of a -1 and J determined by a linear fit of 
m p a in terms of (m T fl) 2 using data at (m n a) 2 pa 
0.3-1.1 (see discussions in the text). Errors are 
statistical only. 
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observed when we introduce the clover term to 
the quark action. For both the P-C and R-C cases, 
the data points drop significantly and lie on top 
of the phenomenological curve within errors. The 
same trend is seen both for rajy jm p and m&/m p . 

A similar effect has been observed in quenched 
QCD by the UKQCD Collaboration in simula- 
tions with the P-C combination at f3 = 5.7-6.2 
(a,- 1 pa l.l-2.5GeV) (§, and also in a study 
with improved gluon actions and the clover and 
D234 quark actions at coarse lattice spacings of 
a,- 1 pa 0.5-0.7GeV @. It is therefore more natu- 
ral to associate the origin of the effect to valence 
quarks rather than to dynamical sea quarks in 
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Figure 2. Static quark potential for the P-W, R-W, P-C, and R-C actions at m^/mp m 0.8. Scales are 
normalized by the lattice spacing determined from m p in the chiral limit. 



full QCD. Nonetheless, changing the gluon and 
quark actions one at a time, we have been able 
to see clearly a decisive role played by the clover 
term in improving the spectrum in a quantitative 
detail in the context of full QCD. In this regard, 
improving the gluon action has much less effect. 

In Table || we compile our results for the J pa- 
rameter, J = my dmv I dm 2 PS at mv/mps = 1.8 
0|. Here again the clover term brings the val- 
ues into better agreement with the experimental 
value of 0.48(2). 

Another interesting feature in our hadron mass 
data is that they exhibit a negative curvature in 
terms of 1/K toward the chiral limit. This is in 
contrast to quenched QCD where hadron masses 
(mass squared for pseudo scalar mesons) are gen- 
erally well described by a linear function of 1/K 
up to quite heavy quark. The curvature is re- 
duced if hadron masses are plotted against raj, 
but still remains at a significant level, especially 
for the R-C combination. This is possibly a full 



QCD effect due to sea quarks which increasingly 
ordered gauge configurations toward the chiral 
limit. For results with clover quark actions the 
trend may be enhanced by the dependence of csw 
on K due to tadpole improvement. 

The curvature causes a practical difficulty in 
the chiral extrapolation of hadron masses since 
our runs have so far been made with at most three 
values of K. For the estimates of a" 1 in Table || 
we employed a linear fit of hadron masses in the 
measured values of mj, excluding for the R-C case 
the point of heaviest quark mass. 

In this connection we note that a comparison 
such as in Fig. [I] and for J in Table || should be 
made at the same lattice spacing in physical units. 
As we see in Table || our estimate of a -1 shows a 
spread of some 20-30% depending on the actions. 
Additional runs are being conducted in an effort 
to match the lattice spacing more precisely. 

Another point to note is that an agreement 
with the phenomenological formula of Ref. p3] 
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Figure 3. Lattice spacing a in GeV 1 as a function of (m^a) 2 , for P-W, R-W, P-C, and R-C actions at 
m^/rrip w 0.8. Note the difference in the scale of horizontal axis for R-C. 



examined in Fig. [I] is not an improvement crite- 
rion that follows from theoretical principles, al- 
beit a reasonable one in view of the success of 
the formula for describing the experimental spec- 
trum. Stability of results toward smaller lattice 
spacing has to be checked, which we hope to pur- 
sue in future simulations. 

4. Static quark potential 

We plot typical results for the static quark po- 
tential in Fig. [2] for which m w /m p 0.8. Con- 
version to physical units is made with the lattice 
scale given in Table ||. Different symbols corre- 
spond to potential data measured in different spa- 
tial directions along the vector given in the figure. 

For the P-W action rotational symmetry is 
badly violated. Improving the quark action (P-C) 
we observe that the potential exhibits a much bet- 
ter rotational symmetry. At present how much 
of the improvement is due to the change of the 
quark action is not clear since the lattice spacing 



for the P-C case (a" 1 « 1.3GeV) is about 20% 
smaller than for the case of P-W (a -1 «s l.IGeV), 
from which one expects a 40% reduction in the 
violation of rotational symmetry. The best im- 
provement is achieved with the R-W and R-C ac- 
tions. It is known that the gauge action R signif- 
icantly improves rotational symmetry already for 
the quenched case |p^| . The improvement natu- 
rally carries over to the present case of full QCD. 

With our present statistics the potential can 
be measured with small errors up to a distance 
of r ss 6a w Ifm. In this distance range we 
do not observe flattening of the potential due 
to pair creation and annihilation effects. In fact 
the potential can be well described by a linear 
plus Coulomb form or + a/r. We extract the 
string tension a by fitting potential data to this 
form. We then use the phenomenological value 
a = (440MeV) 2 to convert results to an estimate 
of a for each value of K. 

In Fig. S we compare the value of a obtained 
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in this way (filled circles) with that from m p = 
770MeV, also calculated for each value of K (open 
circles). For the latter quantity, results obtained 
after an extrapolation of m p a to the chiral limit 
are also shown. This extrapolation is linear for 
a when we take (m^a) 2 as the horizontal axis, as 
adopted in Fig. ||. We observe that the two esti- 
mates converge to a consistent value in the chi- 
ral limit for the P-C and R-C combinations, while 
an apparent deviation of order 40% and 20% are 
indicated for the cases of P-W and R-W respec- 
tively. 

Of course we expect the discrepancy observed 
for the last two cases with the Wilson quark ac- 
tion to disappear in the continuum limit. The 
agreement found for the P-C and R-C cases show 
that the clover term helps improve the consis- 
tency of the two determinations of the physical 
scale of lattice spacing already at a -1 rj 1.3GeV. 

The static quark potential at short distances is 
directly relevant to the spectrum of heavy quark 
systems. A recent NRQCD calculation carried 
out on gauge configurations generated with the 
Kogut-Susskind action for sea quark and the pla- 
quette gluon action at (3 = 5.6 shows a significant 
deviation between the scale determined from the 
heavy and light hadron spectrum fl"i"|| . It would 
be interesting to check for our P-C and R-C cases 
if the lattice scale determined from heavy quark 
systems yield results consistent with those from 
<7 and m p . 

5. Conclusions 

Our comparative study has shown that the 
clover term of Sheikholeslami and Wohlert dras- 
tically improves light hadron spectrum already at 
a -1 w l.l-1.3GeV. While improving the gauge ac- 
tion has much less effect in this regard, we expect 
good rotational symmetry of static quark poten- 
tial, which is achieved with the action R, to be 
important for heavy quark systems. 

Our results lead us to believe that a significant 
step forward towards a realistic simulation of full 
lattice QCD, encompassing both heavy and light 
hadrons, can be achieved with the current gen- 
eration of dedicated parallel computers with the 
application of improved actions. We are particu- 



larly encouraged to pursue the use of the combi- 
nation R-C toward this goal. 
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